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Abstract. We have calculated the distribution of work W done on a 1-d harmonic oscillator 
that is initially in canonical equilibrium at temperature T, then thermally isolated and driven 
by an arbitrary time-dependent cyclic spring constant K{t), and demonstrated that it satisfies 
P{W) = exp(/3W)P(— W), where /3 = l/ksT, in both classical and quantum dynamics. This 
differs from the celebrated Crooks relation of nonequilibrium thermodynamics, since the latter 
relates distributions for forward and backward protocols of driving. We show that it is a special 
case of a symmetry that holds for non-cyclic work processes on the isolated oscillator, and that 
consideration of time reversal invariance shows it to be consistent with the Crooks relation. We 
have verified that the symmetry holds in both classical and quantum treatments of the dynamics, 
but that inherent uncertainty in the latter case leads to greater fluctuations in work performed for 
a given process. 
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1. Introduction 

The Crooks relation UJ states that the outcome of the mechanical processing of a system 
according to a prescribed sequence of actions is related to the outcome of a process consisting 
of the reversed sequence. It is a connection between the probability distributions of the amount 
of work W performed on the system in the course of such forward and backward processes, 
Pp{W) and Pb{W), respectively. A forward process might consist of the movement of a piston 
to compress a gas in a cylinder, while the backward process would be the opposite movement to 
expand the gas [2|. The validity of the Crooks relation requires that the system should start out in 
canonical equilibrium at the same temperature T for both processes. The system might maintain 
contact with a heat bath at that temperature during the processing or it could be isolated. The 
relation reads 

i^(W)=exp(j3(W-AF))PB(-W), (1) 

where j3 = l/ksT and AF is the change in free energy of the system associated with the forward 
process, evaluated for example on the basis of the isothermal change in volume of the expanded 
gas. Equation ([T]) states that the probability that the forward process should require an input of 
work W, and the probability of requiring work — W (in other words receiving work from the 
system) during the backward process, are related to each other, but are not in general equal. 
The relation has been shown to hold for a variety of choices of dynamics, though studies 
reveal that it is important to define carefully what is meant by work, particularly for strong 
coupling between a system and its environment EHH. It implies ||5l the Jarzynski equality 
(exp(— j3 iW — AF))) — 1, which in turn leads to (W) > AF for a forward process starting in 
equilibrium, where the brackets indicate an average taken over the probability distribution of 
work done. This, of course, is a statement of the second law of thermodynamics, and the Crooks 
relation, Jarzynski equality and the associated fluctuation relations have received a great deal of 
attention as a result (see, for example |6l|71[8l|9l[l0l|IT]|T2l[l3[M[l5l[I3[l7l[3[Il). 

It has proved valuable to study the Crooks relation in the context of simple examples 
lfT9l |2] |20l Hi] |22] |23] |24l in order to gain insight into its operation in more complex cases, 
and the 1-d harmonic oscillator has proved to be a popular system. The main purpose of this 
paper is to calculate work distributions for an isolated oscillator using a geometric and pictorial 
approach that has, we believe, some intuitive pedagogical value. In doing so, we expose a 
broader symmetry of the distribution of work for such a forward process, which has its origin 
in the simplicity of the dynamics of the harmonic oscillator, and which demonstrates the rather 
special character of this system with regard to its fluctuation behaviour 

Work is performed by prescribing a time-dependent spring constant K{t) during the process. 
If the spring constant varies cyclically in an interval < f < T such that k{t) = K'(O), then 
AF = 0, and moreover if the process takes place under conditions of thermal isolation, then 
no heat is exchanged during the cycle and W = AE, the change in system energy. The Crooks 
relation reduces to Pp{AE) — exp{l5AE)Ps{—AE): a result verified, for example, by Deffner 
and Lutz 1211 l25l . But our approach to solving the classical evolution demonstrates that a 
relation Pp{AE) = exp{l5LAE)Pp{—AE) also holds for a class of forward processes, with the 
parameter L depending on the nature of the process, and with Pp{AE) taking a specified analytic 
form. It is quite compatible with the Crooks relation for a cyclic forward process and its 
backward counterpart, for which L ~ I and Pp{AE) = Pq{AE), as we shall show. We also 
compare the continuous distribution of work arising from a classical treatment with the discrete 
distribution of work that emerges from the quantum treatment of a cyclic process to show that 
these properties are preserved. Fluctuations in the quantum case are broader, as a result of the 
wider range of possible outcomes made possible by the dynamics. 

In the next section we analyse a general work process performed on a classical harmonic 
oscillator, represented in terms of a matrix operation on a system phasor. We evaluate 
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the probability distribution function (pdf) of system energy change for M > 1 independent 
oscillators. We go on to treat the system quantum mechanically in section [3] particularly to 
contrast the widths of the classical and quantum pdfs for cyclic processes. Our conclusions are 
given in section |4] 

2. Classical treatment 

2.1. System phasor and process matrix 

It is convenient to convert the time-dependent spring constant K{t) into a time-dependent natural 
frequency (o{t), so that the Hamiltonian at time t is 



and the equation of motion is x = —oP'x, where p and x are the momentum and position of the 
oscillator, and m is its mass. The response of the system may be illustrated using pictures of 
phase space orbits. An isolated system with constant frequency COq performs clockwise circular 
orbits in a phase space where momentum is normalised by dividing by m(Bt). The square of the 
radius of the orbit is proportional to the initial energy of the oscillator If the spring constant is 
changed at time f,, altering o and the energy, then the system moves onto an elliptical orbit. A 
process consisting of a sequence of shifts in spring constant without contact with a heat bath can 
therefore be visualised as transitions to, and movement along, a set of elliptical orbits in phase 
space. But if the process is cyclic, characterised by a return to the original spring constant, 
the final orbit will be circular Whether the energy of the system has increased or decreased as 
a result of the process then depends on whether the radius of the final orbit is greater than or 
less than the initial radius, respectively. This is illustrated in Figure [T] for a process consisting 
of a shift down and up in spring constant. The implication is that both upward and downward 
changes in energy can be generated. With the exception of some special cases, where motion 
on the intermediate orbit consists of one complete circuit for example, in which case AZi = 
always, a cyclic process clearly produces a pdf that describes both positive and negative AE. 

Let us consider a general variation in a)(f) in the interval < f < T, with (O ^ (Oq — (i){0) for f < 
and (£) = (On = Co{t) for f > t. In the initial and final situations the motion will be sinusoidal, but 
with differing phases and amplitudes in general. We write x{t) = Aexp(ia)of) + A* exp{—i(Oot) , 
where A is a complex phasor representing the phase and amplitude. The process in the interval 
< f < T will then map an initial phasor A onto a final phasor A' = SA, with a matrix 
representation 



where Ai- and Ai are the real and imaginary parts of the phasor A. We shall call S the process 
matrix. 

The energy of an oscillation with phasor A and angular frequency (O can be written as ImaP'A^A, 
so the change in energy brought about by the process is 



where D = con/cOo and / is the unit 2x2 matrix. In terms of magnitude | A | and phase 
= tan-^(i(0)/a)ox(0)): 



H{t) ^ p^/2m + mCO^{t)x^/2, 



(2) 




(3) 




(4) 




(5) 





Figure 1. The system moves clockwise on a circular orbit in phase space until the spring constant 
changes att = tj and it moves onto an elliptical orbit. The cycle is completed by a change in spring 
constant att = tf to return the system to a circular orbit. Depending on the initial location of the 
system on its orbit, and the duration of the process, the radius of the final orbit could be (a) greater 
than or (b) less than the initial radius, illustrating how the energy change can be both positive or 
negative. 



we find that 

A£=£(Ci+C2COs20+C3sin20), (6) 

where E = 2mft)Q |Ap is the initial energy of the oscillator, and the constants Ci, C2 and C3 can 
be expressed in terms of the elements of S: 

Ci = — {a^+b^ + c^ + d^)-l (7) 



2 

Z)2 



C2 = ) (8) 

C3=D^{ab + cd). (9) 

Let us consider a sequence of step changes in frequency from (0„-i to (0„ at times t„ = n8t, such 
that in the interval tn<t <tn+i the displacement is represented by 

x{t) =A„exp(ia)„(?-?„))+A*exp(-ia^(r-r„)). (10) 
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Imposing continuity requirements for x and i at f = f„, we express A„ in terms of A„_ i : 
An = ^ (^A„_i exp (ia)„_i5f) (^1 + j 

+A:_iexp(-i«„_,5f)(^l-^j^. (11) 

The ife/? matrix S„, representing the transformation from A„_i to A„ through the relation 
A„ ~ SnA„^i, has the following structure: 

where c„_i = cos(a)„_i5f) and i„_i = sin(a)„_i5f), corresponding to a rotation through phase 
angle a)„_ 1 5f , followed by a rescaling of the imaginary part of the phasor. The process matrix S 
may then be constructed from the step matrices as 5 = SnSn-i...Si. Note that detS',, = (B„_i/(B„ 
so that detS' = {(On-i/(On) ■ ■ ■ (00/^1) — (Oo/(On =D^^ =ad — bc and this is unity if the process 
is cyclic with con = (Oo- By introducing small step changes in CO, we can consider processes 
consisting of a continuous variation in spring constant. 



2.2. Construction ofP(AE) 



Consideration of all possible initial conditions establishes the pdf of energy change P{AE) due 
to the process. We assume a canonical distribution over points F in phase space now labelled by 
initial energy E and phase angle 9: 



P(r)dr ^ P{9\E)P(E)dEde = —jis-f^^dEdO, 

2k 



(13) 



2k r^" n 

ile"P^£'(Ci +C2cos20 + C3sin20)d£'d0 
Jo 27r 



and hence the average of is 

(A£> = 

= Ci/j3, (14) 

which, through (|9]l, establishes a connection between (A£) and the elements of the process 
matrix S. 

Consider next the pdf <I> of energy change AE for an oscillator with a given initial energy E. It 
is straightforward to write 



7=1 



de 



dAE 



(15) 



where 9i...9k are the values of 9 that satisfy (|6]). We employ AE from (|6]l in the following form: 

(16) 



AE = E(Ci + {cl + C^y^\os2{9 + 9o)) , 



1 /2 

where cos 200 = C2/ (C| + Cl) ' . Noting that, in non-exceptional cases. A; = 4 we write: 



0(A£,£) = - 



d9 



dtSE 



£ (C| + Cf )'^^sin20 



(17) 
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Figure 2. The pdf of energy change AE for an initial energy E, as specified in )18t . 



where <j) = 9 + Oq, and obtain 



^{AE,E) = TT-i {E^iCj + Cl)-iAE^ECi f) ^'^ 



n [{AE - AE^) {AE+ - AE] 



il/2- 



This is vaHd for AE^ < AE < AE+, where AE± = £ ( Ci ± (C^ + d 



-2^1/2 



(18) 



<I> is zero outside 



this range, is symmetric about AE = EC\, and diverges at AE = AE±, as illustrated in Figure 

m 



Now we include the distribution of initial energies E. We write 



P{AE)oc / (I)(A£',£')exp(-j3£')d£', 
Jo 

which for AE > may be cast more explicitly as 



P{AE)' 

while for AE < we use 

P{AE) ' 



A£Z,+ 



(I>(A£',£')exp(-j3£')d£', 



AEL- 



<P{AE,E)exp{-l5E)dE, 



(19) 



(20) 



(21) 



where L± = yOi ± (C, +C3)'^^j . We have assumed C^+C^ > to ensure that both 
positive and negative AE are generated by the process: as we saw earlier this is almost always the 
case for a cyclic process. The integration limits are best understood by consideration of Figure 
[3j The lower integration limit AEL+ for the AE > case is the energy for which A£'+(£'), the 
upper boundary of the range for which <i>{AE,E) is non-zero, is equal to the given AE (shown 
as the dotted line). Similarly, the lower integration limit AEL^ for the AE < case is the energy 
such that A£'_(£') = -AE. 

It is possible to reduce P{AE) to a simple analytic form for this situation. For positive AE we 
write 

P{AE) oc r ^{AE,E + A£L+)e-/^('^+^^+) dE 
Jo 



3A£Li 



exp(-/3£)d£ 



oc e 



((A£'(L_-L+) -£■)£■) 

-/3Ais(L,+L_)/2^^[^^(^^_^_) 



1/2 



(22) 
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Figure 3. Indication of tlie domain of integration over the pdf ^{AE.E) for the case Ci > 0. The 
function itself is sketched by curves similar to Figure|2]for two values of E. 



where Kq is a modified Bessel function of the second kind. The pdf for AE < may be 
constructed in a similar way: 



P(A£)oc / 4>(A£;£ + A£L_)e-'^('^+^^-)d£ 



oc e 



3A£L_ 



exp(-j3£)d£ 



(-£■(£■ -A£'(L+-L_))) 



1/2 



oc e 



-l3AE{L++L-)/2f 



We immediately notice a symmetry of the pdf of energy change: 

P(A£') =exp(j3LA£')P(-A£'), 



(23) 



(24) 



where L = — (L+ This is reminiscent of the Crooks relation, but the distribution on both 

sides describes a forward process, whereas the Crooks relation concerns forward and reverse 
processes. A consequence is 



(25) 



(exp(-j3LA£')) = / P(A£')exp(-j3LA£')dA£' 

P(-A£')dA£' = 1, 
and by Jensen's inequality we deduce that L{AE) > 0. 

We can express P{AE) in terms of process parameters Ci and D since + C| — Cj 
2Ci + 1 - such that 



-2{Cl + C^y^^ 2{C\+2Ci + l-D 



and 



2^l/2 



Cl C3 



2Ci 



2Ci + 1-D2 



2Ci 



C1-CI-C3 2Ci + 1-D2' 



(26) 



(27) 
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implying that the pdf has a form that depends only on the mean energy change (AE) = Ci /j3 
and the ratio D = (Oq/(On, as long as C 



-c]- 



Cf = 2Ci + l-D^>0. 



For simplicity let us now focus our attention on a cyclic process with (On — dXj and hence D = \, 
L = - (L+ + = 1 and L+ - L_ = (1 + 2/Ci)'/2. The earUer inequality imphes that (AS) > 
and C\ > 0. The pdf takes the form 



P{AE) 



)3exp(j3A£/2) ( P\AE 



PiAE) 



1/2N 



(28) 



where a normalisation constant has been inserted. This form is consistent with (AE) = 
J^^AE P{AE)dAE . It is also consistent with the classical limit of the work distribution of a 
quantum harmonic oscillator for a cyclic process obtained by Deffner and Lutz Ii2ri . 



2.3. Example cases 

The analysis may be illustrated for a cyclic process consisting of a step change in frequency 
from coo to 0)1 at f = 0, followed by another jump from Oi to Oq at f = t. The process matrix S 
takes the form 



5 = 



1 
^ 

(Bo 



Cl -Si 
Si Cl 



1 
^ 

COl 



(29) 



where ci = cos (Oi x and = sin (Oi T which implies that 



Cl = sin ^'^1 '^^ I 2 



(30) 



The results of a Monte Carlo simulation of such a process with Oi = Icoq and (OiT = 7i/2, 
with initial states selected from a canonical distribution with j3 — 1, are shown as a histogram 
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Figure 5. Probability distribution functions of energy cliange AE for ten (solid curve in inset), 
tliree (points), two (dotted) and a single oscillator (dashed), when driven by a work cycle with 
mean energy change {AE) = 2/j} per oscillator 

of energy changes AE (continuous line) in Figure |4] The points are obtained from (l28l l using 
the appropriate value Ci =9/8. The correspondence between the analytical results and the 
simulation is apparent. 

We next calculate Pm{AE), the pdf of energy change forM independent oscillators undergoing 
a given cyclic process, using the iterated convolution operation Pm{AE) = /"„Pi(AZi — 
x)PM-i{x)(ix, where Pi{x) is synonymous with P{x). Since P2{—AE) = f^^Pi{—AE — 
x)Pi{x)dx = /r„^'(-A£'+x)P(-x)djc = /"^/"(AS -x)exp(j3 (-AS +x))P(x)exp(-^x)djc = 
exp(— j3AZi)P2(A£'), the symmetry of the pdf for a single oscillator is retained for a system of 
two oscillators. By iteration, it can be shown that 

/'M(A£)=PM(-A£)exp(j3A£). (31) 

For a cyclic process with Ci = 2 we show numerically generated pdfs in Figure|5]for one, two, 
three and ten oscillators. All of them satisfy (ISTT i in spite of their varying shape as M increases. 

These symmetries of the work distribution, however, are not immediately equivalent to the 
Crooks relation, which refers to forward and backward processes. But it has been remarked 
that the harmonic oscillator is rather special ifTSl |22I and that fluctuation relations of a specific 
kind emerge. Since the Crooks relation states \hatPp{AE) = Pb A/?) exp()3 A/i) for an isolated 
cyclic process it must be the case that Pp{AE) = Ps{AE): the distribution of energy change for 
the oscillator is the same whether we process the system according to a forward cyclic sequence 
(o{t) or the reverse sequence a)(f ) = co{t — t). The special nature of the oscillator allows us to 
understand this in the following way. 

The process matrices for the forward and backward cycles, Sp and 5b respectively, are related 
by SpSS^S = I where 



nr.) 



(32) 
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is a velocity inversion matrix since it transforms phasor A into A*, and these represent oscillator 
configurations with the same position but opposite velocities. The relation is simply the 
statement that a deterministic forward process, followed by velocity inversion, the backward 
process and another velocity inversion, should restore the initial state of the system in phase 
space. Hence, if 



a b 
c d 



then Sn 



d b 
c a 



(33) 



and clearly, both matrices have the same value of C\ according to (|9]l, in spite of corresponding 
to quite different transformations of the initial phasor A. We showed in ( l28T l that the distribution 
of energy change is the same for two cyclic processes with equal values of Ci and (A/i), and 
hence PfIAS) ^PbIA^). 

We noted with reference to Figure |2] that if 



-3 ^ ^1 

or negative values of AZs emerges. For C\ > 0, we find AE. 
take the form 



C? < C? then a distribution with only positive 



> and the distribution would then 



P{AE)' 



A£L_ 



^{AE,E)exp{-l5E)dE, 



(34) 



for positive AE and P(AE) = otherwise, which again can best be understood with reference to 
Figure |3] However, this does not appear to reduce to a simple form in general and a symmetry 



about the point AE = is obviously absent. However, a special case for C| + q = Cf does 



simplify since it corresponds to L 

p{ae: 



oo and L+ = 1 /(2Ci). We find that 
e-P'^dE 



i^EL+ (AE (E - 

PAEL+ / ^ 



AEL 

1/2 



ICinAE 



))l/2 
g-i3A£/(2Ci) 



(35) 



An example process where this applies is a step up from coq to coi > COq: S is simply the right 
hand component matrix in ( |29] | such that 



and the condition 



(36) 




PpiAE) = 



(37) 



for AE > and zero otherwise. For the reverse process consisting of a step down from Oi to (Oq 
the pdf is 

1/2 

expl 



Pb(A£) 



nico 



0)^) AE 



jSco^AE 



©2 - «2 



(38) 



for AE < and zero otherwise. The subscripts identify these as pdfs for a forward process and 
its backward counterpart, and they satisfy 



Pf{AE) _(0o _f Pco^AE 



Pb{-AE) 



— exp 

0)1 



cot 



exp 



PcofAE 



(Of 



mi 



(39) 



= exp(j3(A£-AF)), 

as required by the Crooks relation, where AF = /3^'ln(o)i/o)o) is the free energy change in 
the forward process. They resemble the work distributions for a harmonic oscillator under 
isothermal conditions for this process lITSl . 
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8 10 



Figure 6. Probabilities of energy change AE (in units of /itOo) for a step up-step down cyclic 
work process. Grey bars represent Pq{AE) in j42t . derived from a quantum treatment, and black 
bars correspond to the classical pdf (28), integrated over a range hcog each side of the allowed 
quantised values of AE. The lower plot at Phm = 1 illustrates the similarity of treatments at 
high temperature in contrast to the upper plot at j}tiai = 5 (or low temperature) where the range 
of fluctuations is broader in the quantum treatment compared with the classical. The quantum 
histograms satisfy Pq(A£) = exp(j3A£')Pq(— AE): the classical counterparts do not though they 
are based on underlying continuous pdfs that do. 



3. Quantum treatment 



We expect the pdfs of energy change for a quantum treatment of forward and backward 
processes for an isolated oscillator to satisfy the Crooks relation, as has been demonstrated by 
Deffner et al 11211 125]| . However, our interest in this section is in the symmetry in the distribution 
of energy change for a forward cyclic process. We employ the treatment of a 1-d quantum 
oscillator driven by an arbitrary a){t) provided by Ji et al |26|. According to this approach, the 
familiar eigenfunctions i/„(x) of a 1-d harmonic oscillator with frequency Co{Q) ~ COq evolve 
into 



1 / (0() 



2"n\ \h7tg-{t) 



1/4 



COo 



-dt' 



exp 



(coo- 



i^o(f)),^2 



2hg-it) 



0)0 



lo mg^{t') 

where the Hn are Hermite polynomials and the functions (f ), go(0 ™d g+{t) satisfy 
g-{t)^ -2go{t)/m 
gQ{t) ^mCO^{t)g^{t)-g+{t)/m 
g+{t)^2m(0^{t)go{t), 



(40) 



with initial conditions ^-(0) — 1/m, go{0) —0 and ^+(0) — m(OQ 



(41) 

The pdf of energy change 



Pq{AE) for an arbitrary cyclic process is then straightforward to calculate. We write 



PqijtKOo) = Y.J1 \TUt)\^Pn5k-n.j, 
n=Ok=0 



(42) 
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(a) (6) AE 



Figure?, (ci) The evolution of functions mg_ (?) (thin solid curve), go (?) (dashed) and g+(;)/m£Bj 
(dotted) under time-dependent frequency driving to(f)/£Oo = (l + sin(2cOo;) + sin^(£Oof)) (thick 
solid curve), (b) The probability distribution Pq(A£') of energy change AE (in units of SoJo) for 
jihoo = 1 is plotted in the form f{AE) = exp(— /3A£'/2)Pq(A£') to demonstrate the symmetry 
/(A£) = /(-A£). 



where p„ °^ exp(— (n + ^/2)hcool5) is the initial canonical probability for state n, and 7i-„(f) is 
the transition amplitude, which up to an unimportant phase is given by 




Parity considerations dictate that 7i^„ is zero unless k and n differ by an even number, and hence 
Pq{jh(0(,) is zero unless /' is even. 

We study the step up and down process that was considered classically in section 12.31 namely 
0) = C(Jo for f < and t > T, and co = (Oi = 2c(Jo in the interval <t <T with ©i T = ;r/2. The 
pdfs of energy change under the quantum dynamics for cases where jihcoo is equal to 1 and 5 
are shown in Figure|6] Both pdfs satisfy the relationship Pq(A£') = exp(j3A£')Pq(— AZs). The 
quantum pdfs may be contrasted with the classical counterparts by integrating P{AE) in (l28T l 
over the range (j — l)hcoo < AE < {j + l)hcoo and comparing the result with P(^{jhcoo). We see 
that the classical and quantum treatments coincide rather well for jihcoo = 1, but that differences 
emerge for jihcoo = 5. The latter is evidently a low temperature regime and the energy change 
brought about by the process is distributed more broadly due to the relatively more substantial 
quantum fluctuations. 

The step up-step down process has a very simple time-dependence of Co{t). We 
investigate driving the system with the more complicated frequency history co{t) — 
(Oo (l + sin(2a)of) + sin^(a)of)) over the interval < COQt < n such that the g-functions evolve 
with time as in Figure|3a). The resulting distribution Pq(AZi) for jS/iOJo = 1 is shown in Figure 
ITJ/j) in a form that demonstrates the symmetry Pq(A£') — exp(/3 A£')fq(— A/i). 
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4. Conclusions 

We have calculated the probability distribution function P{AE) of energy change brought about 
by taking a 1-d harmonic oscillator, initially in thermal equilibrium, through a process of 
external work while isolated from the environment. If the process is cyclic, then the pdf 
extends over positive and negative AE and exhibits a symmetry P{AE) = exp(/3A£)P(— AZi) 
with respect to the reversal of the sign of the energy change. This is reminiscent of the 
Crooks relation, which is also satisfied by the system, but is distinct, since it involves a 
forward process only and not its reverse. But the symmetry is also a special case of the 
result P{AE) — e\p{l5LAE)P{—AE) that is valid for a class of non-cyclic processes for this 
system. These symmetries are a consequence of the simple dynamics of an oscillator. We have 
demonstrated that the symmetry is retained if the system consists of M independent oscillators 
subjected to the same process, even though the Pm{AE) take a variety of forms. The symmetry 
is also retained when an oscillator undergoing a cyclic process is treated quantum mechanically. 
At high temperatures the results of the classical and quantum treatments are similar, while at 
low temperatures there is a relative broadening of the pdf in the quantum treatment, as would be 
expected from the inclusion of quantum uncertainty in the dynamics. 

The understanding of nonequilibrium thermodynamic processes has advanced tremendously 
in the last decade or so as a result of the development of fluctuation relations and particular 
identities such as the Crooks relation and the Jarzynski equality fT5ll3l. The explicit calculation 
of probability distribution functions satisfying these relations is a challenging task, but has 
considerable pedagogical value. The harmonic oscillator has been a popular system for such 
activity, but it is a rather special case in that further symmetries emerge that are not present 
in general. Using an approach based on phasors and a geometrical consideration of phase 
space trajectories and their weighting in canonical averages, we have demonstrated some of 
this richness. It offers a contrast to approaches offered elsewhere 11211 l22l l23l l24l |251 with the 
intention that it might make an additional contribution to this understanding. 
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